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ABSTRACT 

'>  We  incorporate  the  effects  of  material  temperature  sensitivity  and  heat 
conduction  into  an  infinite  band  type  analysis  of  shear  localization.  Full 
account  is  taken  of  finite  geometry  changes,  but  inertial  effects  arc  neglected.  An 
energy  balance  is  written  between  homogeneously  deforming  bands  in  a  manner 
that  models  conditions  in  our  recent  finite  element  study  of  shear  band 
development  from  internal  inhomogcneitics  in  a  solid  subject  to  plane  strain 
compression.  The  band  analysis  requires  specification  of  an  imperfection 
amplitude  and  a  length  scale  over  which  heat  conduction  effects  are  significant. 
These  are  chosen  to  match  results  of  a  finite  element  analysis.  The  predictions  of 
the  simple  band  analysis  and  the  results  of  full  finite  element  calculations  arc 
then  compared  for  a  wide  range  of  material  properties  for  both  isotropic  and 
kinematic  hardening  characterizations  of  the  flow  potential  surface.  The  predicted 
dependence  of  the  onset  of  localization  on  material  properties  such  as  strain 
hardening  and  strain  rate  sensitivity  is  the  same  for  both  types  of  analysis.  1'. 
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1.  INTRODUCTION 

Shear  bands  frequently  accompany  the  large  plastic  straining  of  materials. 
The  large  localized  strains  in  a  shear  band  often  precipitate  fracture.  When 
shear  bands  do  not  lead  to  fracture  the  localized  shearing  greatly  affects 
subsequent  plastic  deformation.  Hence,  shear  bands  have  a  dual  significance: 
as  a  precursor  to  fracture  and  as  a  mechanism  for  large  strain  plastic 
deformation.  In  metals  subject  to  high  rates  of  loading,  shear  bands  can  form 
as  the  result  of  a  thcrmomcchanical  instability.  The  same  metal  may  also 
undergo  shear  localization  at  very  low  strain  rates,  where  thermal  effects  are 
negligible,  with  some  other  mechanism  causing  the  localization. 

In  this  paper,  we  present  a  method  for  predicting  the  onset  of  shear 
localization  in  rate  and  temperature  dependent  solids.  Localization  is 
regarded  as  an  instability  which  develops  as  the  result  of  the  growth  of  an 
initial  inhomogcncity  in  a  narrow  plane  (the  imperfection  band)  of  an 
otherwise  homogeneous  material.  We  carry  out  our  analysis  withi  a 
framework  that  regards  localization  as  a  material  instability  (Hill  [1],  Thomas 
[2],  Mandcl  [3],  Rice  [4]),  which  has  been  used  to  analyze  the  influence  of  a 
variety  of  constitutive  features  on  localization,  c.g.  Rice  [4],  Hill  and 
Hutchinson  [5],  and*  Nccdlcman  and  Rice  [6],  The  general  framework  for 
impcrfcction-bascd  localization  studies  for  time  independent  solids  has  been 
presented  by  Rice  [4],  which  follows  the  shear  band  bifurcation  methodology 
inherent  in  Marciniak  and  Kuczynski’s  [7]  analysis  of  localized  necking  in  thin 
sheets.  Studies  of  localization  along  these  lines  have  been  carried  out  by 
Hutchinson  and  Tvcrgaard  [8]  and  Sajc,  Pan  and  Nccdlcman  [9).  The  effects 


m 
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of  material  strain  rate  sensitivity  have  been  included  by  Pan,  Saje  and 
Nccdlcman  [10]. 

The  present  analysis  for  localization  builds  upon  the  existing  framework  for 
rate  dependent  solids  in  [10]  by  incorporating  temperature  sensitivity  into  the 
formulation.  Attention  is  focused  on  the  response  in  the  imperfection  band  and 
in  a  thin  strip  of  material  at  a  characteristic  distance  from  it.  The  conditions 
under  which  a  highly  localized  pattern  of  deformation  in  the  imperfection  band 
can  emerge  relative  to  an  imposed  homogeneous  deformation  in  the  material  strip 
arc  sought. 

The  onset  of  localization  depends  critically  on  the  details  of  the  assumed 
constitutive  law.  It  is  well  known  that  yield  surface  curvature  plays  an 
important  role  in  the  localization  process.  The  classical  clastic-plastic  solid  with  a 
smooth  yield  surface  is  quite  resistant  to  localization,  Rudnicki  and  Rice  [11]. 
However,  for  a  kinematically  hardening  solid  with  a  smooth  yield  surface, 
localization  is  possible  at  achievable  strain  levels  in  the  three-dimensional  shear 
localization  context,  Hutchinson  and  Tvergaard  [8],  This  effect  has  also  been 
observed  in  the  plane  stress  sheet  necking  context,  Tvergaard  [12]  and  Nccdlcman 
and  Tvergaard  [13]. 

Analyses  arc  performed  using  a  viscoplasticity  theory  which,  in  the  rate 
independent  limit,  corresponds  to  J2  flow  theory  with  combined  isotropic  and 
kinematic  hardening:  The  formulation  incorporates  the  effects  of  thermal 
softening,  strain  hardening,  strain  rate  sensitivity  and  heat  conduction.  The 
momentum  and  energy  balance  equations  arc  coupled  via  the  temperature 
dependent  flow  stress  and  the  internal  heat  generation  resulting  from  plastic 
deformation,  and  arc  solved  simultaneously. 

Full  account  is  taken  of  finite  geometry  changes,  but  inertial  effects  arc 
neglected.  Attention  is  confined  to  a  range  of  strain  rates,  from  say  10'^  see'*  to 


103  see"',  where  material  rate  sensitivity  is  the  main  time  effect  and  higher  strain 


consideration. 
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2.  CONDITION  FOR  LOCALIZATION 

An  infinite  planar  solid  which  contains  a  thin  planar  imperfection  band  is 
shown  in  Fig.  1.  The  imperfection  band  is  denoted  by  the  region  labelled  'b\  and  a 
thin  strip  of  the  surrounding  material  located  at  a  distance  1  from  the  imperfection 
band  is  labelled  's'.  The  unit  normal  vector  to  both  the  imperfection  band  and  the 
material  strip,  denoted  by  v,  is  initially  oriented  at  an  angle  <P0  from  the  compression 
axis.  The  width  of  both  regions  is  denoted  by  w,  and  is  assumed  to  be  sufficiently 
small  so  that  the  properties  and  fields  in  them  may  be  considered  to  be  homogeneous. 
Material  properties  and  field  quantities  in  these  regions  arc  identified  with  a 
superscript  or  subscript  'b  '  or  's'. 

The  reference  configuration  for  the  Lagrangian  formulation  corresponds  to 
the  initial  unstressed  state  with  the  material  points  identified  by  the  Cartesian 
coordinates  x1.  In  the  current  deformed  state,  the  material  points  arc  at  x1  +  u1 
where  u‘arc  the  displacement  components,  and  F  =  I  +  3u/3x. 

The  solid  is  deformed  quasi-staticaliy  to  give  the  deformation  rate  Fs  in  the 


material  strip.  As  a  result  of  the  initial  imperfection,  the  current  value  of  the 
field  quantities  inside  the  imperfection  band  will  generally  differ  from  those 
outside  it.  The  conditions  of  compatability,  equilibrium  and  balance  of  energy 
must  be  satisfied. 

The  requirement  of  continuity  of  the  velocity  field  takes  the  form  (Rice 


•b  •  s 

Fjj  =  F  - j  +  qj  vj 


where  (')  denotes  the  material  time  derivative  and  v  is  the  unit  normal  to  the 
imperfection  band  in  the  undeformed  state.  The  tensor  q  represents  the 


current  deformation  non-uniformity  which  has  accumulated  as  a  result  of  the 


initial  imperfection. 
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Thc  condition  of  incremental  equilibrium  between  the  two  regions 
requires  that 

(n[/  -  n‘J)  vj  =  0  (2.2) 

where  n1-1  arc  the  contravariant  components  of  the  nominal  stress  tensor.  It  is 
assumed  that  the  stress  rate  in  (2.2)  is  expressible  in  the  form 

n'j  =  K‘jki  Fik  +  A'J  .  (2.3) 

The  moduli  KljkJ!  and  A‘j  arc  rate  and  temperature  dependent. 

Combining  (2.1)  through  (2.3)  yields 


vivk<1* 


) 


Fikvi  +  (A 


•j 

s 


Ah)  V: 


(2.4) 


which  is  the  evolution  equation  for  q. 

Heat  conduction  effects  arc  modelled  assuming  that  the  time  rate  of 
change  of  energy  in  the  imperfection  band  due  to  conduction  is  proportional 
to  the  temperature  gradient  (T^  -  Ts)/i.  Accordingly,  energy  balance  between 
the  two  regions  results  in  evolution  equations  for  Ts  and  T^,  given  by 


PcpwTs  =  7-(Tb-Ts)  +  Xwts  :  dj  (2.5) 

2k  j,  p 

PCpwTb  -  —  OVTs)  +  *wT  :  db  .  (2.6) 

In  the  equations  above,  p  is  the  mass  density,  cp  is  the  specific  heat  at  constant 
pressure,  k  is  the  thermal  conductivity,  T  is  the  Kirchhoff  stress  tensor,  and  dp  is 
the  plastic  part  of  the  rate  of  deformation  tensor  d,  which  is  the  symmetric  part 
of  F  F"1.  Heat  flowing  out  of  the  imperfection  band  flows  to  either  side  of  it; 
hence  the  factor  of  two  in  (2.6).  The  fraction  of  plastic  stress  work  which  is 


>  .  ►  ^ V  -  .  . 


>  *•  A  A  J** 


converted  to  heat  is  denoted  by  x,  and  is  typically  in  the  range  of  0.85  to  0.95  for 

ductile  metals  (Taylor  and  Quinney  [14]). 

The  heat  transfer  condition  incorporated  in  (2.5)  and  (2.6)  is  taken  to 

model  the  behavior  in  our  recent  finite  element  study  of  shear  localization  [15] 

from  a  doubly  periodic  array  of  internal  inhomogencitics.  In  those 

calculations  heat  conduction  was  accounted  for  within  each  cell  surrounding  an 

inhomogcncity,  but  there  was  no  heat  flux  across  cell  boundaries. 

In  a  planar  analysis,  (2.4)  through  (2.6)  constitute  a  system  of  four 

differential  equations  which  govern  the  state  of  the  system'.  Localization  occurs 

when  one  of  the  components  of  the  rate  of  deformation  tensor  within  the 

imperfection  band  becomes  infinite.  In  the  rate  independent  limit,  A  =  0,  and 

ij  ki 

localization  occurs  when  the  determinant  of  (IC^  v-vj.)  van>shcs  (Rice  [4]). 

Two  dimensionless  groups  characterize  the  thermal  response, 

_  k  *  °o 

t  -  - —  ,  H  (2.7) 

Pc p  wl€n  P‘pT0 

Here,  € n  is  a  reference  strain  rate,  TQ  is  a  reference  temperature,  and  aQ  is  a 
reference  stress.  Uniform  temperature  behavior  results  as  T  approaches 
infinity,  and  adiabatic  conditions  arc  obtained  as  T  approaches  zero.  The 
characteristic  length  (k/pcp  en)  governs  the  extent  over  which  the  effects  of 

lz 

heat  conduction  arc  significant.  We  note  that  the  length  scale  (wi)  entering 
the  definition  of  non-dimensional  parameter  T  is  different  from  the  length 
scale  used  to  define  a  similar  non-dimensional  parameter  in  our  finite  element 
study  of  shear  localization  in  [15],  The  rate  of  internal  heat  dissipation  is 


controlled  by  the  parameter  0. 
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3.  CONSTITUTIVE  RELATIONS 


The  flow  potential  surfaces  are  taken  to  be  concentric  cylinders  centered 
about  a  stress  state  a  in  stress  space.  The  radius  of  the  flow  potential  surface 
is  denoted  by  oF  and  is  given  by 

f  3  — r* 

aF  =  S :  S  j  (3.1) 

where  S  =  S-ct  The  tensor  S  is  the  Kirchhoff  stress  deviator,  given  by 
S  =  T  —  1/3t  :  I,  where  I  is  the  identity  tensor.  Here,  the  kinematic  and 
isotropic  hardening  plastic  flow  rules  are  based  on  Kirchhoff  stress  rather  than  on 
Cauchy  stress.  As  long  as  clastic  strains  remain  small,  there  is  little  difference 
between  the  two  formulations. 

The  rate  of  deformation  tensor  is  expressed  as  the  sum  of  clastic  and 
plastic  parts  by 


d  =  dc  +  d  .  (3.2) 

The  clastic  part  is  given  by 

dc  =  L'1  :  f  (3.3) 

where  ^  is  the  Jaumann  rate  of  Kirchhoff  stress  and  L  is  a  fourth  order 
tensor  representing  the  clastic  stiffness  of  the  material.  The  plastic  part  of 
the  rate  of  deformation  tensor  is  given  by 


dP  = 


37  - 

2^7  s  ‘ 


Here,  £  is  specified  by  the  power  law  relation 


£°  [ [  1  -B(T-T0)]h( e ) ] 


where  € 0  is  a  reference  strain  rate,  TQ  is  a  reference  temperature,  m  is  the  strain 
rate  hardening  exponent  and  6  specifics  the  thermal  softening  characteristics  of 


the  material.  This  assumed  linear  temperature  dependence  is  representative  of 
carbon  steels  in  the  range  0  -  500°C,  c.g.  Stakcr  [16],  although  at  elevated 
temperatures  the  flow  stress  is  generally  a  non-linear  function  of  temperature. 

NVe  employ  combined  isotropic  and  kinematic  hardening,  Gocl  and  Malvern 
[17],  with  h(?)  and  the  evolution  equation  for  a  given  by 

h(I)  =  X  o0[  1  +  F/€0]N  +  (1-X)  a0  (3.6) 

A  2  P 

ft  -  —  (1  -  X)b  d  (3.7) 

m  °o  _  N-l 

[1-0(T-TO)]  N  —  [1  +  €/€0]  (3.8) 

€o 

Here,  oQ  is  a  reference  stress,  €Q  is  a  reference  strain,  N  is  the  strain  hardening 
exponent  and  X  is  a  constant  ranging  between  zero  and  unity.  Purely  isotropic 
hardening  corresponds  to  X  =  1,  while  purely  kinematic  hardening  is  given  by 
X  =  0.  The  evolution  equation  for  b,  (3.8),  is  taken  so  that  at  constant  strain  rate 
and  temperature,  the  response  for  proportional  loading  is  independent  of  whether 
the  change  in  flow  potential  surface  is  described  by  isotropic  hardening  or  by 
kinematic  hardening,  see  Nccdlcman  and  Tvergaard  [13]. 

An  expression  for  the  Jaumann  rate  of  Kirchhoff  stress  is  obtained  in 
terms  of  rate  quantities  by  combining  (3.2)  through  (3.4), 

a  3e 

r  =  L  :  d - L  :  S  (3.9) 

2  °F 

Although  there  is  no  explicit  yielding  in  this  formulation,  for  small  m  there 
is  an  effective  yield  point  in  that  plastic  strain  rates  arc  very  small  when  the 
numerator  is  less  than  the  denominator  in  (3.5). 
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4.  FINITE  ELEMENT  IMPLEMENTATION 


4.1  Rate  Tangent  Modulus 


In  order  to  increase  the  stable  step  size  over  that  of  an  explicit  Euler  time 
integration  scheme,  Peirce  ct  al.  [18]  employ  a  formulation  which  results  in  a 
tangent  modulus  type  method  for  solving  the  governing  rate  equations.  The 
method  is  explicit  in  that  no  iterations  arc  required.  In  this  method  the  effective 
uniaxial  plastic  strain  rate  €  within  an  increment  is  expressed  as  a  linear 
combination  of  its  rates  at  times  t  and  t  +  At  by 


e  =  (1  -  0)€t  +  6et+At  (4.1) 

where  0  is  a  parameter  ranging  from  zero  to  unity.  The  plastic  strain  rate  at  time 


t  +  At  is  approximated  by 


[de  I  3e  3e  I  ■  "I 

=  €f+  —  Op  +  —  £  +  -  T  At  . 

1  L3oFlt  F  3e  t  3Tlt  J 


An  expression  for  ap  is  obtained  by  differentiating  (3.1)  and  using  (3.7)  and  (3.9), 


dF  =  - —  L  :  S  :  d  -  [3G  +  (1  -  \)b]? 

2  Op 


where  G  is  the  clastic  shear  modulus  and  b  is  given  by  (3.8). 
Combining  (4.1)  through  (4.3)  and  solving  for  7  yields 
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Substitution  of  (4.4)  into  the  expression  for  the  Jaumann  rate  of  Kirchhoff 
stress  given  by  (3.9)  results  in 


f  =  C  :  d  -  —  L  :  S 

2of 

(4.8) 

where 

f  3  f  - 

C  =  L  -  1^—  J  q  (L  :  S)  (L  :  S) 

(4.9) 

The  expression  for  £  given  by  (4.8)  may  be  used  to  compute  the  stress  rates  t1J 
needed  for  the  computation  of  the  moduli  K1-*^  in  (2.4).  The  details  of  this 
computation  arc  given  by  Nccdlcman  [19]. 

4.2  Solution  of  the  Coupled  System 

Equations  (2.4)  through  (2.6)  constitute  a  coupled  system  of  four  differential 

equations  which  may  be  solved  to  obtain  the  state  of  the  system  at  time  t  +  At. 

These  equations  may  be  effectively  uncoupled  by  introducing  an  estimate  for  the 

■  c  •  c 

temperature  rates  at  time  t  +  0At,  say  Ts  and  T^,  into  (2.4).  This  is  accomplished 

by  expressing  T  as  a  quadratic  function  of  time  based  on  its  values  at  the 

previous  three  points  in  time  at  which  the  solution  has  been  obtained,  and 

•  e 

extrapolating  to  obtain  T  . 

Following  the  computation  of  the  qj  using  (2.4),  the  deformation  gradient 
-b 

rate  Fjj  is  computed  from  (2.1).  Then  the  deformation  gradient,  stresses  and 
internal  variables  in  the  imperfection  band  and  the  material  strip  are  updated 
using  a  simple  linear  incremental  updating  scheme.  Next,  the  heat  generation  rate 
XT  :  dP  is  computed  in  both  regions,  and  their  temperatures  arc  calculated  from 
(2.5)  and  (2.6).  The  stepsize  is  reduced  if  the  computed  and  the  estimated 
temperature  rates  differ  by  more  than  an  allowable  tolerance. 

When  the  solution  of  the  coupled  system  at  time  t  +  At  has  been  completed, 
this  solution  becomes  the  new  current  state.  The  procedure  is  then  repeated  to 
obtain  the  solution  for  the  next  increment. 


5.  NUMERICAL  RESULTS 

The  material  properties  used  in  the  calculations  arc  chosen  to  be 
representative  of  a  4340  steel  studied  experimentally  by  Hartley  [20],  and  arc 
identical  to  those  used  in  the  full  two-dimensional  finite  element  calculations 
of  shear  banding  by  LeMonds  and  Needlcman  [15].  These  properties  arc 
specified  by  aQ  =  1250  MPa,  E  =  200,000  MPa,  v  =  0.30,  «0  =  0.003,  N  =  0.08, 
m  =  0.01,  e0  =  0.001  see'1,  p  =  7833  kg/m3,  cp  =  465  J/kg°C,  k  =  54  \V/m°C 

(in  the  analyses  involving  heat  conduction),  0  =  0.0016  per  °C,  and  T0  =  20 °C. 

These  values  comprise  a  reference  set  of  material  properties  which  arc  used 
for  comparison  of  critical  strain  predictions  of  the  flow  localization  analysis  to 
the  critical  strains  obtained  from  the  finite  element  calculations  in  [15]. 
Non-dimensional  parameters  which  depend  on  aQ  arc  o0/E  =  0.00625  and 
n  =  1.1714.  The  initial  undeformed  state  is  stress  free  and  at  a  uniform 

temperature  specified  by  TQ. 

The  imperfection  is  prescribed  by  specifying  a  value  for  the  initial  flow 
stress  in  the  imperfection  band  which  is  slightly  different  from  that  in  the 
surrounding  material.  We  define  a  measure  of  the  amplitude  of  the 

imperfection  with  the  parameter  f  =  Og/c^.  Values  of  f  less  than  unity  arc 

considered,  implying  that  the  material  in  the  imperfection  band  is  initially 
softer  than  the  surrounding  material. 

The  solid  is  compressed  under  conditions  of  plane  strain  at  a  constant 
nominal  strain  rate  4n,  with  the  compression  axis  parallel  to  the  x3  axis.  This 
deformation  is  obtained  by  prescribing  F22  =  €nF22,  Fj2  =  F|t  =  0,  and  hjX  =  0. 
The  stress  rate  boundary  condition  may  be  used  in  conjunction  with  (2.3)  to 


The  sensitivity  of  the  critical  strain  to  the  imperfection  strength  f  is  shown 
in  Fig.  2  for  an  isotropically  hardening  material  deforming  under  adiabatic 
conditions  with  €n/«0  =  5x10^  (en  =  500  see**);  therefore  T  =  0  in  (2.7)  and 
X  =  1  in  (3.6)  and  (3.7).  The  parameter  ecr  denotes  the  critical  strain,  which  is 
the  lowest  value  of  the  localization  strain  for  all  possible  initial  orientations  of 
the  imperfection  band.  The  localization  strain  is  the  value  of  the  maximum 
principle  logarithmic  strain  in  the  material  strip  at  the  point  when  7^/7.  =  10. 
The  critical  strain  is  very  sensitive  to  the  imperfection  strength. 

In  subsequent  calculations,  we  use  a  value  of  f  =  0.9994  which  yields  the 
critical  strain  ecr  =  0.150  observed  in  the  finite  element  calculation  of  the 
isotropically  hardening  solid  deforming  under  adiabatic  conditions  in  [15]. 
Although  the  identification  of  the  initiation  of  localization  in  the  finite 
element  calculations  is  somewhat  arbitrary,  the  strain  accumulations  in  the 
band  become  very  rapid  at  the  indicated  critical  strain.  The  initial  orientation 
of  the  band  is  <P0  -  51.7°,  and  its  orientation  at  localization,  given  by 
tan  ^cr  =  (Fi1tan0o/F22)cr,  is  <l>cr  =  43.1°.  This  result  agrees  well  with 
<Pcr  =  42.8°  observed  in  the  finite  element  calculation. 

The  effect  of  the  thermal  parameter  T  on  the  localization  strain  is  shown  in 
Fig.  3  for  a  kinematically  hardening  solid  deformed  at  en/«0  =  5x10^.  The 
adiabatic  and  uniform  temperature  limits  can  be  specified  to  within  five  percent 
error  (relative  to  the  numerical  results)  by  T  <  0.06  and  T  >  175,  respectively. 
The  critical  strains  at  the  adiabatic  and  uniform  temperature  limits  arc 
€ cr  =  0.170  and  0.873,  respectively,  with  corresponding  initial  band  orientations 
of  <1>Q  =  52.4°  and  76.0°. 

Finite  element  calculations  in  [15]  for  a  kinematically  hardening  solid  with 
heat  conduction  effects  were  carried  out  for  €n/e0  =  5x10^,  5x10“*  and  5x10^, 
which  correspond  to  €n  =  500  see'*,  50  see'*  and  5  see**.  In  each  ease, 


localization  occurred  rapidly  once  a  shear  band  initiated.  The  parameter  T  in 
(2.7)  characterizes  the  thermal  response  of  the  solid,  with  T  =  0  specifying  the 
adiabatic  limit  and  T  **  "  specifying  the  uniform  temperature  limit.  The  value 
of  T  which  results  in  localization  at  the  critical  strains  at  which  the  shear  bands 
formed  in  the  finite  element  calculations  in  [15]  is  shown  in  Table  1  for  each  of 
the  normalized  strain  rates  considered.  These  values  of  T  will  be  used  in 
subsequent  calculations  which  include  the  effects  of  heat  conduction. 

Several  finite  element  calculations  using  the  numerical  method  in  [15]  have 
been  carried  out  here  at  €n/«0  =  5x10^  (en  =  500  scc’^)  for  kinematically 
hardening  solids  with  material  properties  which  arc  significantly  different  from 
those  used  in  the  previous  calculations.  The  critical  strains  predicted  by  the  flow 
localization  analysis  agree  quite  well  with  the  strains  at  which  shear  bands 
formed  in  the  finite  element  calculations,  as  shown  in  Table  2. 

Figures  4  through  7  depict  €cr  as  a  function  of  the  material  properties  6,  N, 
m  and  n,  respectively,  for  a  kinematically  hardening  solid  deforming  at 
€n/€Q  =  5x10-*  with  X  =  3.37  in  (2.7).  In  these  figures,  the  results  of  the  finite 
element  calculations  reported  in  Table  2  arc  plotted  with  the  symbol  ‘X’.  The 
agreement  between  the  predictions  of  the  flow  localization  method  (using  a 
constant  imperfection  amplitude  f  =  0.9994)  and  the  results  of  the  finite  element 
analyses  arc  quite  good  over  a  range  of  considerable  variation  in  the  material 
properties. 

The  localization  strain  predicted  in  the  temperature  independent  limit  (0  -  0 
in  (3.5))  in  Fig.  4  is  finite,  although  this  is  not  readily  apparent  from  the  figure. 
The  localization  strain  in  this  ease  is  very  high  due  to  the  weak  imperfection 
(f  =  0.9994).  A  stronger  imperfection  will  permit  localization  at  a  strain  which  is 
physically  plausible.  For  example,  ecr  =  0.50  when  f  =  0.991.  If  an  isotropic 
hardening  surface  is  used,  then  the  material  is  so  resistant  to  localization  when 
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8  =  0  that  an  imperfection  specified  by  f  =  0.878  is  necessary  in  order  to  obtain 
£cr  =  ®-50. 

The  critical  strain  shown  as  a  function  of  the  strain  hardening  exponent  N 
(see  (3.6))  in  Fig.  5  exhibits  a  maximum  value  emax  =  0.296  at  N  =  0.22.  The 
material  stiffness  and  the  rate  of  internal  heat  generation  both  increase  with 
increasing  N.  Strain  hardening  tends  to  stabilize  the  deformation,  while  softening 
due  to  increased  heating  is  destabilizing,  Clifton  [21].  These  competing  effects 
account  for  the  observed  maximum  critical  strain. 

Fig.  8  shows  the  dependence  of  the  localization  strain  on  the  imperfection 
band  orientation  for  several  values  of  the  thermal  softening  parameter  8  in  (3.5) 
for  the  ease  of  kinematic  hardening  with  X  =  3.37  in  (2.7).  The  solid  curves 
correspond  to  the  initial  band  orientation,  and  the  dashed  curves  correspond  to  the 
angle  at  localization.  The  imperfection  band  orientation  at  localization  always 
falls  below  77/4,  which  is  consistent  with  the  predictions  of  bifurcation  analyses 
for  plane  strain  compression,  Hutchinson  and  Tvergaard  [8]. 

The  jump  in  hydrostatic  stress  across  the  band  interface  is  shown  in 
Fig.  9  for  several  values  of  the  strain  hardening  exponent  N  in  (3.7)  for  a 
kinematically  hardening  solid  deformed  at  en/€0  =  5*10-*  with  X  =  3.37  in 
(2.7).  For  N  >  0.17,  the  hydrostatic  stress  in  the  imperfection  band  is  less  than 
that  in  the  surrounding  material  throughout  the  deformation  history,  whereas 
for  more  lightly  hardening  materials  the  jump  in  hydrostatic  stress  is  negative 
initially  and  then  increases.  For  the  ease  where  N  =  0.05,  the  maximum 
principal  logarithmic  strain  is  0.224  inside  the  band  and  0.167  outside  the  band 
when  the  increase  in  the  jump  in  hydrostatic  stress  takes  place.  This 
suggests  the  possibility  that  at  large  accumulated  shear,  the  hydrostatic  stress 
in  the  band  may  become  positive  even  though  the  nominal  stress  state  is  one 
of  plane  strain  compression.  This  is  of  interest  in  regard  to  ductile  failure 


mechanisms  since  void  growth  is  driven  by  a  positive  hydrostatic  tension  (see 
c.g.  Rice  and  Tracey  [22],  McClintock  [23]). 

The  effect  of  the  mixed  hardening  parameter  X  in  (3.6)  and  (3.7)  on  the 
localization  strain  is  shown  in  Fig.  10  for  a  solid  compressed  at  €n/€0  =  5x10^ 
(«n  =  5  see"')  with  T  =  25.34  in  (2.7).  Localization  occurs  at  €cr  =  0.554  with 
tp0  =  69.3°  for  kinematic  hardening,  and  at  £cr  =  0.93  with  4>Q  =  80.3°  for 
isotropic  hardening.  The  effective  stress  in  shear  is  lower  for  the  kinematic 
hardening  surface,  and  the  higher  curvature  of  that  surface  relative  to  the 
isotropic  hardening  surface  results  in  a  greater  increment  in  shear  strain  for  a 
given  increment  in  stress.  In  the  finite  element  calculation  in  [15]  with  isotropic 
hardening  at  this  imposed  strain  rate,  the  solid  was  deformed  until  e  =  0.90 
without  any  indication  of  shear  band  formation. 
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6.  CONCLUSIONS 

The  results  of  the  present  analysis  indicate  that  shear  localizaiton  in  rate  and 
temperature  dependent  solids  is  strongly  influenced  by  a  variety  of  factors.  Which 
factors  arc  most  influential  depends  on  the  complete  description  of  the  material  -  its 
properties  and  constitutive  behavior  -  and  the  imposed  deformation  rate.  For  the 
kinematically  hardening  solid  specified  by  our  reference  set  of  material  parameters, 
and  deformed  in  plane  strain  compression  at  €n  =  500  see*1,  a  five-fold  increase  in 
the  localization  strain  occurs  between  the  adiabatic  and  uniform  temperature  limits, 
identifying  heat  conduction  as  the  primary  stabilizing  effect.  The  main  destabilizing 
factor  in  this  ease  is  thermally-induced  softening:  a  physically  plausible  value  of  the 
thermal  softening  parameter  (0  =  0.0016  per  °C)  permits  a  relatively  weak 
imperfection  (f  =  0.9994)  to  result  in  localization  at  ecr  =  0.217.  The  results  also 
illustrate  the  effect  of  the  interaction  of  the  material  properties.  For  example, 
whether  strain  hardening  is  stabilizing  or  destabilizing  depends  on  the  degree  of 
thermal  softening.  Furthermore,  as  strain  rates  decrease,  the  effects  of  thermal 
softening  diminish,  and  the  description  of  the  flow  potential  surface  plays  a  greater 
role  in  the  localization  process.  In  the  analyses  with  en  =  500  see'*,  the  particular 
description  of  the  combination  of  isotropic  and  kinematic  hardening  is  only  a 
secondary  effect,  altering  the  localization  strain  by  no  more  than  twenty  percent. 

Using  physically  plausible  values  for  the  imperfection  amplitude  f  and  the 
thermal  parameter  T,’thc  flow  localization  method  yields  critical  strains  which  arc 
in  very  close  agreement  with  the  results  of  the  full,  two-dimensional  finite  element 
calculations  of  shear  banding  by  LeMonds  and  Nccdlcman  [15].  Our  present 
calculations  illustrate  that  the  use  of  constant  values  of  f  and  T  result  in 
reasonably  accurate  critical  strain  predictions  over  a  fairly  broad  range  of 


material  properties. 


-17- 


ACKNOWLEDGEMENTS 

The  support  of  this  work  through  ARO  Grant  DAAG29-85-K-0003  is 
gratefully  acknowledged.  The  computations  reported  on  here  were  carried  out  at 
the  Brown  University,  Division  of  Engineering,  Computational  Mechanics  Facility. 
The  acquisition  of  this  facility  was  made  possible  by  grants  from  the  U.S. 
National  Science  Foundation  (Grant  ENG78-19378),  the  General  Electric 
Foundation,  the  Ford  Motor  Company  and  the  Digital  Equipment  Corporation. 


REFERENCES 


[1]  Hill,  R.,  “Acceleration  Waves  in  Solids,”  J.  Mech.  Phys.  Solids,  Vol.  13,  p. 
10-16,  1962. 

[2]  Thomas,  T.Y.,  Plastic  Flow  and  Fracture  in  Solids.  Academic  Press,  New 
York,  1961. 

[3]  Mandcl,  J.,  Rheology  and  Soil  Mechanics,  cds.  J.  Kravtchenko  and  P.M. 
Siricys,  Springcr-Vcrlag,  1966. 

[4]  Rice,  J.R.,  “The  Localization  of  Plastic  Deformation,”  in  Theoretical  and 
Applied  Mechanics,  cd.  W.T.  Koitcr,  North-Holland,  Amsterdam,  pp.  207-220, 
1976. 

[5]  Hill,  R.  and  Hutchinson,  J.W.,  “Bifurcation  Phenomena  in  the  Plane  Tension 
Test,”  J.  Meek.  Phys.  Solids,  Vol.  23,  pp.  239-264,  1975. 

[6]  Nccdlcman,  A.  and  Rice,  J.R.,  “Limits  to  Ductility  Set  by  Plastic  Flow 
Localization,”  in  Mechanics  of  Sheet  Metal  Forming,  eds.  D.P.  Koistinen  and 
N.M.  Wang,  Plenum  Press,  New  York,  1978. 

[7]  Marciniak,  Z.  and  Kuezynski,  K.,  “Limit  Strains  in  the  Processes  of 
Stretch-Forming  Sheet  Metal,”  hit.  J.  Mech.  Sci„  Vol.  9,  1967. 

[8]  Hutchinson,  J.W.  and  Tvergaard,  V.,  “Shear  Band  Formation  in  Plane 
Strain,”  Int.  J.  Solids  Struct.,  Vol.  17,  pp.  451-470,  1981. 

[9]  Sajc  M.,  Pan,  J.  and  Nccdlcman  A.,  “Void  Nuclcation  Effects  on  Shear 
Localization  in  Porous  Plastic  Solids,”  Int.  J.  Fracture,  Vol.  19,  pp.  163-182, 
1982. 


[10]  Pan,  J.,  Sajc  M.  and  Nccdlcman  A.,  “Localization  of  Deformation  in  Rate 
Sensitive  Porous  Plastic  Solids,”  hit.  J.  Fracture,  Vol.  21,  pp.  261-278,  1983. 

[11]  Rudnicki,  J.W.  and  Rice,  J.R.,  “Conditions  for  the  Localization  of 
Deformation  in  •Pressure-Sensitive  Dilatant  Materials,”  J.  Mech.  Phys.  Solids, 
Vol.  23,  pp.  371-394,  1970. 

[12]  Tvergaard,  V.,  “Effect  of  Kinematic  Hardening  on  Localized  Necking  in 
Biaxially  Stretched  Sheets,”  Int.  J.  Mech.  Sci.,  Vol.  20,  pp.  651-658,  1978. 

[13]  Nccdlcman,  A.  and  Tvergaard,  V.,  “Limits  to  Formability  in  Rate  Sensitive 
Metal  Sheets,”  in  Mechanical  Behavior  of  Materials,  cds.  J.  Carlsson  and  N.G. 
Ohlson,  Pergammon,  pp.  51-65,  1984. 

[14]  Taylor,  G.I.  and  Quinney,  H.,  “The  Latent  Energy  Remaining  in  a  Metal 
After  Cold  Working,”  Proc.  Roy.  Soc.  London  143,  1934. 


j.  v  /  tl 


[15]  LcMonds  J.  and  Ncedlcman,  A.,  “Shear  Localization  in  Rate  and 
Temperature  Sensitive  Materials,”  to  be  published. 

[16]  Stakcr,  M.R.,  “The  Relation  Between  Adiabatic  Shear  Instability  Strain  and 
Material  Properties,”  Acta  Metall.,  Vol.  32,  pp.  683-689,  1981. 

[17]  Gocl,  R.P.  and  Malvern,  L.E.,  “Biaxial  Plastic  Simple  Waves  With  Combined 
Kinematic  and  Isotropic  Hardening,”  J.  Appl.  Mech.,  Vol.  37,  1970. 

[18]  Peirce,  D.,  Shih,  C.F.,  and  Ncedlcman,  A.,  “A  Tangent  Modulus  Method  for 
Rate  Dependent  Solids,”  Computers  and  Structures,  Vol.  18,  pp.  875-887, 
1984. 

[19]  Ncedlcman,  A.,  “Finite  Elements  for  Finite  Strain  Plasticity  Problems,”  in 
Plasticity  of  Metals  at  Finite  Strain:  Theory.  Computation,  and  Experiment, 
edited  by  E.H.  Lee  and  R.H.  Mallctt,  pp.  387-436,  Stanford  University,  1982. 

[20]  Hartley,  K.A.,  “Temperature  Measurements  During  the  Formation  of  Shear 
Bands  at  High  Rates  of  Deformation,”  Ph.D.  Thesis,  Brown  University, 
1986. 

[21]  Clifton,  R.J.,  "Adiabatic  Shear  Banding,”  Chapter  8  in  Materials  Response  to 
Ultra  High  Loading  Rates.  National  Materials  Advisory  Committee, 
NMAB-356,  1980. 

[22]  Rice,  J.R.  and  Tracey,  D.M.,  "On  the  Ductile  Enlargement  of  Voids  in 
Triaxial  Stress  Fields,”  J.  Mech.  Phys.  Solids,  Vol.  17,  pp.  201-217,  1969. 

[23]  McClintock,  F.A.,  "A  Criterion  for  Ductile  Fracture  by  the  Growth  of  Holes," 
J.  Appl.  Mech.,  Vol.  35,  pp.  363-371,  1968. 


TABLE  1 


Values  of  the  non-dimensional  thermal  parameter  T  in  (2.7)  which  result  in 
localization  at  the  critical  strain  which  marks  the  point  of  shear  band  initiation 
in  the  finite  element  calculations  with  kinematic  hardening  presented  in  [15],  The 
adiabatic  limit  is  specified  by  T  "*  0  and  the  uniform  temperature  limit  is  given 
by  (  -  The  imperfection  amplitude  is  f  =  0.9994.  The  values  of  all  other 
material  properties  arc  specified  in  the  text. 


TABLE  2 


Critical  strain  predictions  of  the  flow  localization  method  compared  with  the 


critical  strains  at  which 

shear  bands 

were  observed 

initiated  in  finite  element 

calculations 

carried  out  here  using  the 

analytical 

method  presented  in 

[15].  The 

imperfection 

amplitude  is 

f  =  0.9994. 

The  values 

of 

all  other  material 

properties 

arc  specified 

in  the  text. 

Parameter 

Value 

T 

\ 

ccr-FEA 

£cr-FLOC 

0 

0.00105 

3.37 

0 

0.330 

0.347 

6 

0.0020 

3.37 

0 

0.168 

0.168 

N 

0.20 

3.37 

0 

0.295 

0.296 

N 

0.30 

3.37 

0 

0.277 

0.286 

m 

0.025 

0 

I 

0.193 

0.200 

m 

0.002 

3.37 

0 

0.175 

0.150 

m 

0.025 

3.37 

0 

0.283 

0.293 

n 

0.75 

3.37 

0 

0.345 

0.359 

n 

1.50 

3.37 

0 

0.174 

0.164 
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Fig.  1 


Fig.  2 


Fig.  3 


Fig.  4 


Fig.  5 


The  imperfection  band  (b)  and  strip  of  surrounding  material  (s) 

separated  by  the  characteristic  length  1.  Plane  strain  compression 

occurs  in  the  x  direction. 

Localization  strain  €cr  vs.  the  imperfection  amplitude  f  =  Oq/o^  for 
an  isotropically  hardening  solid  deformed  at  €n/£0  =  ^x10^ 
(«n  =  500  see"*).  The  value  of  the  Fourier  modulus  in  (2.7)  is  (,  =  0, 
which  corresponds  to  the  adiabatic  limit.  The  values  of  all  other 
material  properties  arc  specified  in  the  text. 

Localization  strain  ecr  vs.  the  non-dimensional  parameter  T  in  (2.7) 
for  a  kinematically  hardening  solid  deformed  at  en/€0  =  5x]0"* 

(en  •=  500  see'1),  with  the  imperfection  amplitude  f  =  0.9994.  The 
value:  of  all  other  material  properties  arc  specified  in  the  text. 

Localization  strain  €cr  vs.  the  thermal  softening  parameter  13  in  (3.5) 
for  a  kinematically  hardening  solid  deformed  at  4n/e0  =  5*10^ 

(«n  =  500  see'1),  with  the  imperfection  amplitude  f  =  0.9994  and 
T  =  3.37  in  (2.7).  The  values  of  all  other  material  properties  arc 

specified  in  the  text.  Points  labelled  with  an  ’X1  correspond  to 

results  of  finite  element  calculations. 

Localization  strain  ecr  vs.  the  strairt  hardening  exponent  N  in  (3.6)  for 
a  kinematically  hardening  solid  deformed  at  en/£0  =  5x10"* 
(en  =  500  see'1),  with  the  imperfection  amplitude  f  =  0.9994  and 
T  =  3.37  in  (2.7).  The  values  of  all  other  material  properties  arc 

specified  in  the  text.  Points  labelled  with  an  ’X’  correspond  to 

results  of  finite  element  calculations. 


Localization  strain  €cr  vs.  the  strain  rate  sensitivity  exponent  m  in 
(3.5)  for  a  kinematically  hardening  solid  deformed  at  «n/e0  =  5x10-* 
(«n  -  500  sec*'),  with  the  imperfection  amplitude  f  =  0.9994  and 
X  =  3.37  in  (2.7).  The  values  of  all  other  material  properties  are 
specified  in  the  text.  Points  labelled  with  an  'X'  correspond  to 
results  of  finite  element  calculations. 

Localization  strain  «cr  vs.  the  non-dimensional  parameter  n  in 
(2.7),  which  measures  the  magnitude  of  the  rate  of  internal  heat 
generation,  for  a  kinematically  hardening  solid  deformed  at 
4n/€0  =  5x10-*  (4n  =  500  sec*1),  with  the  imperfection  amplitude 
f  =  0.9994  and  T  =  3.37  in  (2.7).  The  values  of  all  other 
material  properties  arc  specified  in  the  text.  Points  labelled  with 
an  'X'  correspond  to  results  of  finite  clement  calculations. 

Localization  strain  €  vs.  the  imperfection  plane  orientation  angle  0 

( -  initial  angle,  -  -  -  angle  at  localization)  for  a  kinematically 

hardening  solid  deformed  at  «n/e0  =  5x10^  (tn  =  500  sec"1),  with 
the  imperfection  amplitude  f  =  0.9994  and  T  =  3.37  in  (2.7).  The 
values  of  all  other  material  properties  arc  specified  in  the  text. 

Jump  in  hydrostatic  stress  across  the  imperfection  band  interface, 
normalized  with  respect  to  aQ  in  (3.6),  vs.  e,  the  maximum  principle 
logarithmic  strain  outside  the  band.  The  solid  hardens  kinematically 
and  is  deformed  at  €n/«0  =  5x10^  («n  =  500  see*'),  with  the 

imperfection  amplitude  f  =  0.9994  and  X  =  3.37  in  (2.7).  The 
values  of  all  other  material  properties  arc  specified  in  the  text. 

Localization  strain  €cr  vs.  the  mixed  hardening  parameter  X  in  (3.6) 
and  (3.7)  for  a  solid  deformed  at  <n/«0  =  5x10-*  («n  =  5  sec"1), 
with  the  imperfection  amplitude  f  =  0.9994  and  X  =  25.34  in  (2.7). 
The  values  of  all  other  material  properties  arc  specified  in  the  text. 


ro  0.01 


0  01 

~®  N  =  0.05 

t 0=5Z.3°  \ 

•  Qn-o.io 

’/'0  =  56.2o  1 

0  02 

~©  N  =  0.I5 

'I'q  '58.1°  1 

'  ®N  =  0.25 

f  o  *58.3° 

)  03  ■ 

— 1 — 1 — i — i — 1 — i —  i 

-i-  J  [  r  ,  | 

0 

0  05 

1  1  1  1  }  | 

°1  0.15 

e 

FIGURE  9 

